Elementary formula for the Hall conductivity of interacting systems 
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A formula for the Hall conductivity of interacting electrons is given under the assumption that the 
ground state manifold is iV gs -fold degenerate and discrete translation symmetry is neither explicitly 
nor spontaneously broken. 



I. 



INTRODUCTION AND RESULTS 



The Hall conductivity in the integer (IQHE) and frac- 
tional (FQHE) quantum Hall effect |jj is observed to be 
quantized. This quantization can be explained by ex- 
pressing the Hall conductivity c H in terms of a topo- 
logical invariant C that takes integer values through the 
relation 



= —uC 



(1.1) 



where v is the filling fraction of the single-particle Lan- 
dau levels For example, the relation between the 
quantized number C and the non- degenerate many-body 
ground state wave functions \^>{4>,p)) obeying twisted 
boundary conditions parametrized by the pair of angles 
< </>, p < 27T for a gas of electrons confined in two- 
dimensional (position) space and subjected to a uniform 
magnetic field is 
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(1.2) 

when the cyclotron energy dominates over the character- 
istic electron-electron interactions while the characteris- 
tic spatial variations of any one-body potential felt by 
the electrons are much longer than the magnetic length 

SI- . 

In this paper, we derive a formula for the Hall conduc- 
tivity that applies to electrons moving in d-dimensional 
space as long as (i) there exists a gap between the N gs - 
fold degenerate ground state manifold and (ii) discrete 
translation invariance holds 0. We show that the Hall 
conductivity ave raged over the degenerate ground states 
is given by Eq. (5.11), which reduces to 



d 2 k F(k) ri(fc) 



(1.3) 



in two-dimensional space if restricted to the case when 
the many-body ground states are exclusively built out of 
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a single Bloch band with the Berry curvature F(k). All 
the many-body corre latio ns in Eq. ([h|) are encoded by 
ii(fe) defined in Eq. (5.8a), i.e., the expectation value of 



the occupation number operator of the Bloch momentum 
k from the Brillouin zone fi averaged over the -/V gs -fold 
degenerate many-body ground states. 

Equation (1.3) reproduces the following results. 

Cases of the IQHE or the Chern band insulator. When 
N gs = 1 and the many-body ground state is the Slater de- 
terminant made of all available Bloch states in the band, 
fi(fe) = 1 for all k G il and 



a H = e - J d 2 k F(k) =:jC 



(1.4) 



with C the Chern number of the band (in particular, 
the Berry curvature is the constant F(k) = 1/51 and the 
Chern number is C = 1 for the lowest Landau level). 

Case of the FQHE. When iV gs > 1 is the degeneracy 
expected from a hierarchical ground state in the FQHE 
at the partial filling v of the lowest Landau level, the 
Berry curvature is the constant F(k) = C/Cl with C = 1 
the Chern number of the lowest Landau level and 



e 



d 2 k 

n 



n(fc) 



(1.5) 



Case of the anomalous Hall effect. When N gs = 1 
and the many-bo dy g round state is a Fermi liquid, we 
can interpret Eq. (1.3) as the anomalous Hall conductiv- 
ity assuming the order of limits by which the gap has 
been taken to zero after all other limits have been taken. 
The anomalous Hall conductivity has been derived in the 
noninteracting limit pTo| , but its derivation allows for 
interactions in this paper. 

Equation ( jl.3] ) predicts the following results for frac- 
tional Chern insulators |TT -21 . 



1. The integral over the Brillouin zone of F(k) x fi(fe) 
equals a rational number p/q, since Laughlin's 
gauge argument for quantization Q then applies. 

2. The integral over the Brillouin zone of F(k) x n(fc) 
obeys a sum rule, for one cannot change the rational 
number p/q continuously. 



2 



3. The Hall conductivity er H does not need to be equal 
to the filling fraction v obtained by integrating 
ii(fe)/f2 over whenever the Berry curvature F(k) 
is not uniform over the Brillouin zone. In this case, 
if F-„ and F mav are the minimum and maximum 

III 111 III dX 

of F(k) over the Brillouin zone, respectively, then 

2 2 

j v (0 F min ) <a n < e -v{n F max ) . (1.6) 

The remaining of the paper is devoted to deriving these 
results. 



III. 



LINEAR RESPONSE 



II. 



DEFINITIONS 



We consider N e electrons confined to a box of volume 
V in position space R d where d — 1, 2, • • ■ whose dynam- 
ics is governed by the conserved many-body Hamiltonian 
H Q acting on the many-body Hilbert space $ N . Bound- 
ary conditions have been imposed such that H has a 
countable number of eigenstates. The spectral decompo- 
sition of H is written 

H a = E gs J2\n)(n\+ ]T E m \m)(m\. (2.1a) 

n=l m=7V gs + l 

The ground state eigenenergy E gs is assumed AT -fold 
degenerate. We reserve the integer- valued index n = 
1, ■ • • , N gs for its N gs linearly independent ground states. 
A gap separating E gs from the eigenenergies E m of all ex- 
cited states is assumed, whereby we reserve the integer- 
valued index m = N gs + 1, • • • for the excited states. 
Eigenstates of H are normalized according to 

„_„ (ro|n>=0. (2.1b) 



(n\n') 



(m\m') 



This means that the scalar product of any two states 
from the Hilbert space $ N is dimcnsionlcss. Finally, the 
total electric charge is assumed conserved by H and this 
symmetry is not spontaneously broken. 

Because of the boundary conditions, the set of center 
of mass momenta is countable, i.e., we can resolve the 
identity according to 

1 = E P Q' PqPq'=Sq, Q/ Pq, (2.2) 
Q 

where Q £ M. d denotes a center of mass momentum and 
Pq denotes the projector onto states with this center of 
mass momentum. 

Finally, we shall assume that neither does H break 
discrete translation invariance [Q, i.e., 

[Hv,Pq}=0 (2-3) 

for any center of mass momenta Q, nor is discrete trans- 
lation invariance spontaneously broken. 

We work in units where the electric charge e, the 
Planck constant h, and the speed of light c are all unity. 
Therefore, only dimensions of energy appear in the the- 
ory. 



We are after the linear response of one ground state 
\n) of the many-body Hamiltonian H to adiabatically 
switching on a spatially homogeneous and static electric 
field £. The coupling to £ is described by the Hamilto- 
nian 

H — Ho + H^t), H 1 (t):=-X-£e r > t , (3.1) 

where X = (X i ) is the many-body position operator in 
the Schrodinger picture |^2|, p3| , i = 1, • • ■ ,d labels the 
spatial coordinates, and r\ is a small positive number that 
implements the adiabatic turn-on of £. 

We would like to make explicit the existence of a 
conserved electronic current that couples to the applied 
static and uniform electric field £. To this end, observe 
that we can always write 



£ = -(d t A)(t), A(t):=-t£. 



(3.2a) 



Insertion of Eq. (3.2a) into Eq. (|3.l|) gives 



H = H Q — J ■ A(t) e"' +d t [X- A(t) e r "] (3.2b) 



with the current |2J] 

J :=i[H Q ,X] = d t X. 



(3.2c) 



In effect, we have performed the time-dependent gauge 
transformation 



H(t) :=H -J-A(t) e"* 
H (t) :=e +^(*)<"' 

J(t) :=i\H (t),X 



Hoe -iX-A(t), 



(3.3) 



In anticipation of linear response theory, we are going 
to use the interaction picture that endows any operator 
O in the Schrodinger picture with the time dependence 



O 1 ^) : = e +iH ^0 e 



iH„ t 



(3.4) 



We do the same with the instantaneous ground states 
\n) t and the instantaneous excited states \fh) t of H (t). 
In linear response theory, we approximate the time evo- 
lution of any state from $ N in the interaction picture by 
linearizing in the perturbation 

H\ (t) = J 1 ■ A(t) e" * + d t [X 1 ■ A{t) e" *] . (3.5) 

For example, for one of the degenerate ground states, 



In(t)) 1 := (l-i J dt'Hl{t' 




(3.6) 



Hence, if we rule out any level crossing [see Eq. (3.12) 
and Ref. 



J, then 

t (n\J(t)\n) t = l {n(t)\j l (t)\n(t)) 1 



(3.7) 
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to linear order in the time-dependent perturb ation and 
after making use of the identity [following Eq. ( |3.2cj) 1 



(njj^^ln) = (n\J\ri) = 



(3.8) 



that disposes of the lowest order contribution to the ex- 
pansion. With some intermediary steps involving the 

integral representation J 1 ^) = J dt 2 J l (t 2 ) and 

1 — oo 

the use of partial integration, we arrive for any given 
i = 1, ■•• ,d at the leading order estimate (summation 
convention over the repeated index j = 1, • • ■ , d) 



t( n \Ji(t) \ n )t = i£j 



t t 
dt-L 



dt 2 (n\[jj(t),jj(t 2 )]\n). 

(3.9) 



Finally, insertion of the resolution of the identity in $ N 

i = 5>>M + £l m >H (3.io) 

n m 

delivers 

i~\ 7 u\ \~\ c V- H J i \m){m\J j \n) - (io j) 
t {n\J l (i) \n) t = - i£j 



-EL 



Here, we used that 



(n\J\ri) 



(3.11) 
(3.12) 



for any pair n, n' — 1, • • ■ ,^V gs i n view of Eq. (3.2c). 



Equation (3.12) prevents level crossing among the degen- 



erate ground states. For any i,j = l,---,d, we introduce 



the conductivity tensor crj^ through 
t(n\Ji (t) \n) t 



2nVa^£ y 



(3.13) 



For any pair i ^ j = 1, • • ■ , d, the volume V on the right- 
hand side guarantees that the Hall conductivity tensor 



(n|JJm)(m|J» - 



2tt^ ^ 



■EL, 



(3.14) 



is (i) intensive despite the fact that the left-hand side 
scales with t he sys tem size and (ii) dimensionless if d = 
2. Equation (3.14) is the conventional representation of 
the Hall conductivity tensor. The right-hand side is an 
infinite series which is presumed convergent because of 
the energy denominators. 

For our purpose, it is, however, more useful to trade 
the current for the position operator in the matrix ele- 
ments of th e right-hand side of Eq. (3.14) with the help 
of Eq. ( |3.2c ) and, in turn, dispose of the energy denomi- 
nator 



^ 2ttV 



^£[{n|x ( MH*>>-(**+j)]- 

(3.15) 



However, we must be very careful with the interpretation 
of the matrix element (m\Xj\n). In the thermodynamic 
limit N e ,V — > oo holding the electronic density N c /V 
fixed, the position operator X - is unbounded, its expec- 
tation value in any momentum cigenstate is ill-defined, 
and so is its trace over the Hilbcrt space. Hence, the rep- 
resentation (3.15) of the Hall conductivity tensor is done 
in terms of the difference of two series, each of which is 
divergent in the thermodynamic limit N c , V — > oo hold- 
ing the electronic density N e /V fixed. Our goal is to 
properly regularize the subtraction of two infinities on 
the right-hand side by means of formal algebraic manip- 
ulations. 



IV. 



ALGEBRAIC REGULARIZATION 



For any i ^ j = 1, ■ ■ • , d, instead of the Hall conduc- 
tivity of a single degenerate ground state, we define the 
average 



(n) 



(4.1) 



Because of the normalization (2.1b), we can introduce 
the projector on the ground states 



P g :=J2\n)(n\ 

n 

and the projector on the excited states 



Evidently, the normalization (2.1b) implies that 

P 2 = P g = P+, P 2 = P = PL 
P P = P P = 0, P 



e ' 

P„ = l. 



(4.2a) 



(4.2b) 



(4.2c) 



In terms of the projectors (4.2a) and (4.2b), Eq. (4.1) 
becomes 



Tr [P^X^X^-ii^j)} (4.3) 



where Tr denotes the trace over the Hilbcrt space $ N . 

The discrete translational invariance of the Hamilto- 
nian H guarantees that it commutes with Pq for any 
center of mass momentum Q, Pq H = H Pq. Corre- 
spondingly, any eigenstate of H can be labeled by a 
center of mass momentum and 



p q p s = p s p q> 



Pq P c — P c Pq j 



VQ. (4.4) 



Since we ruled out spontaneous symmetry breaking 
of discrete translation invariance by assumption, any 
ground state must have a vanishing center of mass mo- 
mentum. If so, 



p qP s 



p s Pq 



0. 



VQ^O, 



(4.5) 
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so that, combined with the application 



P Q = p q (p s 



PA = P,+PA P. 



of the resolution of the identity, we deduce that 



P 



Q 



P Q P C 



P, P 



VQ ^ 0. 



(4.6) 



(4.7) 



We now conjure the following additive decomposition 
of the position operator X = (X i ) H , 



X = T + A. 



(4.8) 



On the one hand, the component T i of T = (T { ) shifts 
momentum by an infinitesimal amount in the i = 1, • • • , <i 
direction. On the other hand, the operator A = (AJ 
commutes with any Pq. It is of crucial importance to 
note that the decomposition ( |4.8| ) is not unique ||^, |23| . 

We are now in position to do the following manipula- 
tions for any « = !,-•■ , d. First, we do the decomposition 



P e X i P g = P C (T. +A i )P g 

i<3 T i _ 1 

i<7 



P r lim 



P„ + P„ A, P„. 



(4.9) 



Second, we use the orthogonality P c P g = to dispose of 
the term i/q that would blow up in the limit of q — > 0, 



P. X P. 



lim P r 
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P, + P r A, P„ 



(4.10) 



Third, wc make use of the r esol ution of the identity and 
the o rtho gonality from Eq. ( |2.2| ) together with Eqs. (4.5) 
and (4.7) to infer that 



P, X, P„ 



lim P, 



i<7 



P„ 



Pg pQ=0- 



(4.11) 

The first term on the right-hand side connects two sec- 
tors of the Hilbert space $ N with well-defined center 
of mass momenta differing by the momentum Q = qe i . 
The second term on the right-hand side annihilates any 
many-body state with nonvanishing Q. Henceforth, the 
product 



p„ x. Pi x, p„ = p„ x 4 p, x, p„ 



(4.12) 



on the right-hand side of Eq. (|4.3|) becomes 



P g X { P c X jPg = P g A A j P g - P g A P g A, P g (4.13) 

if i 7^ j (so that P ^ P ) and where we have as- 

sumed that we can freely interchange the limit with the 
evaluation of the products. 

We now insert Eq. ( [4 . 1 3| ) into Eq. (p~3|), 



Tr 



2tt V 

(P g A Pg)(P g A j P g 



P g [A,^.]p g 

(P g A J .P g )(P g AP g 



(4.14) 



The full trace Tr over the Hilbert space $ N is thus re- 
duced to a trace over the ground state manifold. This is 
a finite sum since we have assumed that the ground state 
manifold is a finite dimensional vector space. To dispose 
of the second contribution, we make use of the cyclicity 
of the trace restricted to the ground state manifold. We 
are then left with the finite sum 



2nVN gs 



Tr [PJA^PJ 



(4.15) 



2irV N, 



Equation (4.15) is the main result of this paper. It is 
an algebraic counterpart to the many-body presentation 
of the Hall conductance in terms of a many-body Berry 
phase defined by twisting boundary conditions 



V. BLOCH REPRESENTATION 

To proceed, we need to choose a basis of the Hilbert 
space $ N . We choose the basis that follows from using 
the Fock space ^ spanned by the creation and annihi- 
lation operators xl(k) and XaC 2 ), respectively, whereby 
any such pair labeled by the band index a = 1, • • • , N 
and the Bloch momentum k corresponds to a Bloch state 
\x^ a H^))- We choose the normalization conventions 



{X a (k),xi(k')} = ns(k-k')6 ati 



nvs 



k,k' ^a,b 



(5.1) 



given the volume V in position space set by the infrared 
cutoff (the linear size L 3> o say) and the volume f2 in 
momentum space set by the ultraviolet cutoff (the lattice 
spacing a <C L say). In other words, we have assumed 
that H obeys the additive decomposition 



H, 



Int 



(5.2a) 



Here 



rrBlo 



N 

E 

a=l 



d d k 



e a (k)xi(k) Xa (k) (5.2b) 



where, for any band index a = 1, • • • , N and Bloch mo- 
mentum k, the single-particle eigenvalue s a (k) and the 
single-particle Bloch state |x (^)) are the solution to 
the eigenvalue problem 



K°(k)\x^(k))=e a (k)\ X ^(k)}. 



(5.2c) 



The interacting term -ffp 11 ' is of higher order than two in 
the number of creation and annihilation operators. 

The decomposition of the position operator on the 
right-hand side of Eq. (4.8) can now be understood as 
follows. Assume that the N x N Hermitcan matrix 
%Q lo (fc) has been specified in the basis \tp^(k)) where 
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a = 1, ■ ■ • ,N. We shall call this basis the orbital ba- 
sis. Diagonalization of %Q lo (fc) delivers the Bloch basis 
Ix^C 8 )) where a = 1, •• • ,N. The unitary transforma- 
tion that brings the orbital to the Bloch basis has the N 2 
matrix elements ita (fe) where a, a = 1, • • • , N. Any pair 
of columns or rows from this matrix must be orthogonal, 

ul a >(k) uW(fc) = 5 atbt ui a Hk) u[ a) *(fc) = S a>b . 

(5.3) 

Here and in what follows, the summation convention 
over repeated band a = 1, • • ■ ,N or orbital index a = 
1, ■ • • ,N is implied. For any spatial coordinate i = 
1, • • • , d, any pair a, b = 1, • • • ,N of bands, and any 
Bloch momentum k we define the non-Abelian Berry con- 
nection 



I 'ft (b) \ 

A? (fe) := -iu^*(k) -jjjL- (fe). 



(5.4) 



Hence, for any pair i j = 1, • ■ • , d, we conclude with 
the desired representation 



d d k 



'ftAf\ 



n a6 (fc) (5.11) 



of the Hall conductivity tensor averaged over the degen- 
erate ground states. 



A. Noninteracting band insulator 

For a noninteracting band insulator with the lowest 
N < N bands filled, we have N„ = 1 and 



n ab (fc) = nf (fe) 



0. 



1 < a,b < N, 
otherwise. 



(5.12) 



The operators on the right-hand side of Eq. (4.8) are then 
represented by [p3l 



and 



A = J ^xi(k)A« b (k) Xb (k) 



(5.5) 



(5.6) 



Q 

Consequently, for any i ^ j = 1, • • ■ ,d, 



[A,^] = / ^xl(fe)[A(fe),^(fe)] oi, X6(fe)- (5-7) 



Equation ( ^>.7[ ) suggests that we define the A^ 2 dimen- 
sionless intensive numbers 



n ab (fc) := 



1 



nvN a 



^2(n\xi(k)x b {k)\n). 



(5.8a) 



S s n=l 

With the normalization ( |5.1| ), one verifies that 

< n QQ (fc) < 1 (5.8b) 

for any band a and any momentum k. Insertion of 
Eq. (Oh into Eq. <Kuh yields 



(5.9) 



Due to th e presence of the gap, the Bloch Hamilto- 
nian ( 5.2c ) can be adiabatically deformed to 



H° [ °(k) — H-2P(fc), 



(5.13) 



where -P(fc) is the projector on the single-particl e sta tes 
in the lower bands. The right-hand side of Eq. (5.13) is 
invariant under any unitary transformation of the filled 
bands, an element of the unitary group U(N) which is a 
subgroup of the group of unitary transformati ons U (N) 
that mixes all N bands. The Hall conductivity ( 5.11 ) can 
be written in a form for which this symmetry manifests 
itself as a local U(N) gauge invariancc, 



d d k 
~2T 



(5.14) 



Here, the non-Abelian Berry curvature of the lower bands 
is given by 

FfP{k):=d ki Af{k) + iAf{k)Af{k) -(i++j) (5.15) 

and the indices with " are only summed over the lower 
bands. Note that 

Af(k) Af(k) -(i*+j) = Af(k) Af{k) 

(5.16) 

vanishes if both a and c arc summed over. This allows us 



to recast the expression for the Hall conductivity (5.11) 
in terms of the manifestly U(N) gauge invariant Berry 
curvature (5.14). 



It remains to evaluate with the help of Eq. (^J) the 
commutator [A i (k) 1 A^(k)] ah of the non-Abelian Berry 
connection. It is 

[^(fc), J 4 J (fc)]» 6 =i l-g- Uk)-(i«j). (5.10) 



B. Noninteracting Fermi sea 



Let us now consider a Fermi sea ground state for which 
-ZV gs = 1. Even though our derivation relies on the exis- 
tence of a gap from the outset, we can think of letting 
this gap go to zero at the end, as long as the ground state 
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remains unique. We shall work with a single partially oc- 
cupied band 5=1 for simplicity. With FS C f2 denoting 
the Fermi sea 

f 1 U p T?Q 

n ab (k) ee n ab (fc) = S a , x S h i x i ' ' (5.17) 

V ' y ' 0,1 I 0, G !J \ FS. ; 



The Hall conductivity ( 5.11 ) becomes 



' 2^ 

FSCO 



(5.18) 



This result agrees with the zero temperature limit of the 
Hall conductivity for the anomalous Hall effect || [To| . 



C. Interacting partially filled single Bloch band 

Finally, we assume that the ground states \n) have non- 
vanishing amplitudes only with those Slater determinants 
that are made of single-particle Bl och states from the 
band 5 = 1. The Hall conductivity (3.11) becomes 



d d k 
~2T 



Ofc)n(fc), 



(5.19) 



we conclude that 



^ x inf F£(k) < a i3 < ^ x sup Fg(k) (5.21) 

27T fcco 2-7T 



feen 



for any pair i,j = 1, • ■ ■ , d EfJ . 

When d = 2, we deduce from Laughlin's flux insertion 
argument H that 



_P_ 



(5.22) 



where N gs is the topological ground state degeneracy on 
the 2-torus and p is any integer that does not need to 
be co -prim e with N gs . If we combine Eq. ( 5.22Q with 



Eq. ( 5.21 ), we conclude that 



2tt 



supfy(fc) - inf i^(fc) 
ken ken 



1 

> - x 



N, 



(5.23) 



is a necessary condition for the Hall conductivity to de- 
viate from 



C\ 



^F${k), (5.24) 



where fi(fe) is the occupation number of the state k in 
that single band 5 = 1 that contributes to any one of 
the N gs ground states \n). The right-hand side is invari- 
ant under any U(l) gauge transformation of the U(l) 
Berry connection A 1 (k) that defines the Berry curvature 
(fc). 

Equation (5T9J) sets bounds to the Hall conductivity 
that can arise in an interacting system from a partially 
occupied single Bloch band. If we define the filling factor 



F 11 



d d k 
~0~ 



n(fc), 



(5.20) 



in two dimensions. Such a deviation was discussed for the 
FQHE in a periodic potential |27j and could in principle 
also appear in two-dimensional fractional Chern insula- 
tors Mi. 
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